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We perform a detailed phase-space analysis of various phantom cosmological models, where the 
dark energy sector interacts with the dark matter one. We examine whether there exist late-time 
scaling attractors, corresponding to an accelerating universe and possessing dark energy and dark 
matter densities of the same order. We find that all the examined models, although accepting stable 
late-time accelerated solutions, cannot alleviate the coincidence problem, unless one imposes a form 
of fine-tuning in the model parameters. It seems that interacting phantom cosmology cannot fulfill 
the basic requirement that led to its construction. 



I. INTRODUCTION 



Recent cosmological observations support that the universe is experiencing an accelerated expansion, and that the 
transition from the deceleration phase to the accelerated phase happened in the recent pastil L In order to explain 
this remarkable phenomena, one can modify the theory of gravity, such as /(i?) gravity [2, |3[, Dvali-Gabadadze- 
Porrati model [3| , and string inspired models Q . An alternative direction is to introduce the concept of dark energy 
which provides the acceleration mechanism. The simplest candidate of dark energy which is consistent with current 
observations is the cosmological constant. Due to the lack of a good explanation of the small value of the cosmological 
constant, many dynamical dark energy models were explored, such as a canonical scalar field (quintessence) model 
@, a phantom model that has equation of state parameter w < — 1 7, 8, 9], or the combination of quintessence and 
phantom in a unified model named quintom jicj ]. Moreover, many theoretical studies are devoted to shed light on 
dark energy within some quantum gravitational principles, such as the so-called "holographic dark energy" proposal 

The dynamical nature of dark energy introduces a new cosmological problem, namely why are the densities of 
vacuum energy and dark matter nearly equal today although they scale independently during the expansion history. 
To solve the problem, one would require that the matter density and dark energy density always approach their current 
values independent of the initial conditions. The elaboration of this 'coincidence' problem led to the consideration 
of generalized versions of the aforementioned models with the inclusion of a coupling between dark energy and dark 
matter. Thus, various forms of "interacting" dark energy models [HI, [l3| have been constructed in order to fulfil the 
observational requirements. 

One decisive test for dark energy models is the investigation of their phase-space analysis. In particular, to examine 
whether they possess attractor solutions corresponding to accelerating universes and ratio fidark energy/ ^dark matter of 
the order 1. If these conditions are fulfilled, then the universe will result to that solution at late times, independently of 
the initial conditions, and the basic observational requirements will be satisfied. Although non-interacting quintessence 
[3 (3 j non- interacting phantom [l6| and non-interacting quintom models [l7| admit late-time accelerated attractors, 
they possess fidark energy = 1 and thus they are unable to solve the coincidence problem. 

When dark energy - dark matter interaction is switched on in quintessence, then one can find accelerated attractors 
which moreover give fidark encrgy/^dark matter ~ O(l) fl8l Il9l. I20I [2l| . but paying the price of introducing new problems 
such as .justifying a non-trivial, almost tuned, sequence of cosmological epochs [22j . In cases of interacting phantom 
models [231 ] . the existing literature remains in some special coupling forms with mainly numerical results, which 
suggest that the coincidence problem might be alleviated [ID, . 
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In the present work we are interested in performing a detailed phase-space analysis of the interacting phantom 
paradigm, considering new and general forms of the interaction term Q. In addition, we go be yon d the simplified 
non-local forms of the literature, considering also a local Q-form proportional to a density [2ll. l25l. I26L \2l\ . Doing 
so we find that although late-time accelerated attractors do exist in all the models, as it is expected for phantom 
cosmology, almost all of them correspond to complete dark energy domination and thus they are unable to solve the 
coincidence problem. Only for a narrow area of the parameter space of one particular model, can the corresponding 
solution alleviate the aforementioned problem. 

The plan of the work is as follows: In section[TT]we construct the interacting phantom cosmological scenario and we 
present the formalism for its transformation into an autonomous dynamical system, suitable for a phase-space stability 
analysis. In section ITTT1 we perform the phase-space analysis for four different models, using various interaction terms, 
and in scction lTVl we discuss the corresponding cosmological implications. Finally, section|V]is devoted to the summary 
of the obtained results. 



II. PHANTOM COSMOLOGY 



Let us construct the interacting phantom cosmological paradigm. Throughout the work we consider a flat Robertson- 
Walker metric: 

ds 2 = dt 2 - a 2 (t)dx 2 , (1) 

with a the scale factor. 

The evolution equations for the phantom and dark matter density (considered as dust for simplicity) are: 

p m + 3Hp m = -Q (2) 
p + 3ff(p0 +p lt> ) = Q, (3) 

with Q a general interaction term and H the Hubble parameter. Therefore, Q > corresponds to energy transfer from 
dark matter to dark energy, while Q < corresponds to dark energy transformation to dark matter. For simplicity, 
we take the usual phantom energy density and pressure: 

94, = -\fi + V{(j>) (4) 

= -\fi-V{<j>)t (5) 

where V(<p) is the phantom potential. Equivalently, the phantom evolution equation can be written as: 

Z mutt A 9V M Q *n 

(j) + ZH(j) — — = r. 6 

Op (j) 

Finally, the system of equation closes by considering the Friedmann equations: 

K 2 

H 2 = ^(Pt+Pm), (7) 
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,- 2 

H = — 



y(p0+P0 + P™), (8) 



where we have set k 2 = 8irG. Although we could straightforwardly include baryonic matter and radiation in the 
model, for simplicity reasons we neglect them. 

In phantom cosmological models, the dark energy is attributed to the phantom field, and its equation of state is 
given by 

uu = — . (9) 

Alternatively one could construct the equivalent uncoupled model described by: 

p m + 3H(1 + W m<e ff)p m = (10) 
/30 + 3H(1 + w^ e ff)p<t> — 0, (11) 
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where 

Wm.eff = 1T§— (12) 
W^ff = W ^-J§^~- ( 13 ) 

However, it is more convenient to introduce the "total" energy density ptot = Pm + p<j>, obtaining: 

ptot + 3H(l+w tot )p tot = 0, (14) 

with 

wtot = — ~r — = WthQth, (15) 

P4, + Pm 

where = = ^dark energy Obviously, since p to t — 3H 2 /k 2 , (TT4")) leads to a scale factor evolution of the form 
a(t) oc £ 2 /( 3 ( 1 +'" , *°*)) ; in the constant Wtot case. However, in the late-time stationary solutions that we are studying in 
the present work, Wtot has reached to a constant value and thus the above behavior is valid. Therefore, we conclude 
that in such stationary solutions the condition for acceleration is just w to t < —1/3. 

In order to perform the phase-space and stability analysis of the phantom model at hand, we have to transform 
the aforementioned dynamical system into its autonomous form [LH . Il5j . This will be achieved by introducing the 
auxiliary variables: 

H0 



(16) 



together with M = In a. Thus, it is easy to see that for every quantity F we acquire F = H-^jj. Using these variables 
we obtain: 

f^^W+y 2 , (17) 

— — 

W<j, = 2~- 2' ( 18 ) 

and 

wtot = -x 2 ~ y 2 . (19) 

We mention that relations (fl8|) and (Tl9|) are always valid, that is independently of the specific state of the system (they 
are valid in the whole phase-space and not only at the critical points) . Finally, note that in the case of complete dark 
energy domination, i.e p m — » and f2^ — > 1, we acquire w to t ~ W4, < —1, as expected to happen in phantom-dominated 
cosmology. 

The next step is the introduction of a specific ansatz for the interaction term Q. In this case the equations of 
motion ©, (O and © can be transformed to an autonomous system containing the variables x and y and their 
derivatives with respect to M = In a. The consideration of various Q-aiisatzes i s performed in the next section. 

A final assumption must be made in order to handle the potential derivative that is present in ([6]). The usual 
assumption in the literature is to assume an exponential potential of the form 

V = V exp(-K\(t>), (20) 

since exponential potentials are known to be significant in various cosmological models [lil [l5[ . Note that equivalently, 

1 

' WW) ' 



but more generally, we could consider potentials satisfying A — — d « const (for example this relation is 
valid for arbitrary but nearly flat potentials 28]). 

Having transformed the cosmological system into its autonomous form: 



X' = f(X), (21) 
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where X is the column vector constituted by the auxiliary variables, f(X) the corresponding column vector of the 
autonomous equations, and prime denotes derivative with respect to M — In a, we extract its critical points X c 
satisfying X' = 0. Then, in order to determine the stability properties of these critical points, we expand ([21) around 
X c , setting X = X c + U with U the perturbations of the variables considered as a column vector. Thus, for each 
critical point we expand the equations for the perturbations up to the first order as: 

U' = S • U, (22) 

where the matrix S contains the coefficients of the perturbation equations. Thus, for each critical point, the eigenvalues 
of S determine its type and stability. 



III. PHASE-SPACE ANALYSIS 



In the previous section we constructed the interacting phantom cosmological model, with an arbitrary interaction 
term Q, and we presented the formalism for its transformation into an autonomous dynamical system, suitable for 
a phase-space stability analysis. In this section we introduce four specific forms for Q and we perform a complete 
phase-space analysis. 



A. Interacting model 1 



This specific interacting model is characterized by a coupling of the form [2Cj, |2l|, |29( : 

Q = aHp m . 



(23) 



Thus, inserting the auxiliary variables (|16|) into the equations of motion 
in the following autonomous system: 



([7]) and ((8]), and using |23|) we result 



x = — 3x + 2 X (1 ~ x2 



V 



a 

2x, 



V 



y(l -x 2 - y 2 ) 



Xxy. 



(24) 



The critical points (x c , y c ) of the autonomous system (fM]) are obtained by setting the left hand sides of the equations 
to zero. The real and physically meaningful (i.e corresponding to y > and < f2^ < 1) of them are presented in 
table HI In the same table we present the necessary conditions for their existence. 
The 2x2 matrix S of the linearized perturbation equations writes: 



-1 



(i- 



(>.r : 



■(1 



vl) 



lie 



-V6A - 3x c + ^ 



-,/|A - 3x r 



~3y 2 c - 



Xx r 



Therefore, for each critical point of table HI we examine the sign of the real part of the eigenvalues of H, which 
determines the type and stability of this specific critical point. In table H] we present the results of the stability 
analysis. In addition, for each critical point we calculate the values of w to t (given by relation (TTT)]0 . and of fi^ (given 
by |17p). Thus, knowing Wtot we can express the acceleration condition wtot < — 1/3 in terms of the model parameters. 

In order to present the results more transparently, in fig. [T] we depict the range of the 2D parameter space (a, A) 
that corresponds to an existing and stable critical point B of table HI i.e to negative real parts of the corresponding 
eigenvalues. Furthermore, we numerically evolve the autonomous system (|24l) for the parameters a = 0.5, A = 1.0 
and a = —3.1, A = 0.1, and the results are shown in figure [2l Depending on which region of the parameter-space do 
the chosen parameter-pair belong (see table (p}), the system lies in the basin of attraction of either the critical point 
A or B, and thus it is attracted by the one or the other. This is apparent in figure^ where we present two different 
parameter choices, leading respectively to point B and A. 



B. Interacting model 2 



In this model we consider an interaction term of the form I2C 



„, „2n tt3— 2n n 

a K n p 



(25) 
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Cr. P. 






Existence 


Stable for 


!^ 


W t ot 


Acceleration 


A 


A 

V6 


>A + S 


all a, A 


-a < A 2 + 3 


1 


x 3 


all a, A 


B 


3+a 
V6\ 


, / a (3+ap 
V 3 6A^ 


A 2 > (3+a)2 for > a > 3 
(3 + a) > A 2 > (3+Q)2 for a < 3 


a < -3 
(3 + a) > A 2 > (3+Q)2 


(3+a) 2 a 
3A a 3 


a 
3 


a < -1 



TABLE I: The real and physically meaningful critical points of interacting model 1 and their behavior. 




-60 -50 -40 -30 -20 -10 



a 



FIG. 1: The range of the 2D parameter space (a, A) that corresponds to an existing and stable critical point B of tabled that 
is in the case of interacting model 1. 




x x 



FIG. 2: Phase-space trajectories for interacting model 1. The left panel corresponds to a = —3.1 and A = 0.1, and thus to the 
stable fixed point B of table|I]with (x c , 2/ c )=(-0.41, 0.93). The right panel corresponds to a = 0.5 and A = 1.0 and thus to the 
stable fixed point A of table|T]with (x c , y c )=(-QAl, 1.08). 



This ansatz corresponds to a generalization of interacting model 1. In the following, we take n = 2 for simplicity. 
Inserting P5|) into the cosmological equations @, ©, and ([5]), using the auxiliary variables, we acquire the 
following autonomous system: 

, V6. 2 3 2 a 3 (l+x 2 ~y 2 ) 2 

x = — Ax \y H — x[l — x — y a 

2 y 2 v y 1 2 u x 

y> = -^-\xy + ^y(l-x 2 -y 2 ). (26) 



6 



The real and physically meaningful critical points (x c , y c ) of the autonomous system (|26|) are: 

(Xci = ± 



-OLQ 

1 + a 



A 



Vd 



0), (x c2 = --j=, y C 2 



x 2 

1 + i>- 



where 



and 



(xc3, Uc3 = yl-^- V6\x c3 /3), (x c4 , y c4 = \fl - x 2 c4 - V6Xx c4 /3), 
K - 1)A + VA 2 (a -l) 2 -12a 



(27) 



a; C 3 



X c i 



2\/6a 



-(ao ~ 1)A - v/A 2 (a -l) 2 -12a 



2\/6a 

These critical points arc presented in table |TTJ The matrix H of the linearized perturbation equations is 

-a/6 A - 3x c + eao^i^ 



(28) 



(29) 



. Il±^±a _ 2qo(1 + X 2 _ y2) + f Il±d^ 



- A /|A-3x c 



2 



|Ax c 



Examining the eigenvalues of the matrix H for each critical point, we determine its stability behavior. Finally, in 
table HT1 we also present Wtot-, &<f> and the acceleration condition. Note that the values of Wtot and Q,^ for the critical 
points C and D can be straightforwardly calculated using P7|) and (fS5|) . (j2"9")h but we do not present them since these 
points are unstable and thus non-relevant for our discussion. 



Cr. P. 



Existence 



Stable for 



W tot 



Acceleration 



A 



±J- 



l + QQ 



a = 



Unstable 



No 



B 



A 



1 + 



A- 



all A, Qo 



all A, ao 



3 



all A, ao 



C 



Xc3 



Vcli 



A<-,/ 3(1 + ° n) ,-l<o <0 



Unstable 



+ 2 = 



3(1 + Q ) 



/311+aolf 
-a 



A < 



"0 



1/2 < a < 

1 < a < -1/2 



D 



a; C 4 



J/c4 



Unstable 



-^04 + y c 



■>>- 



V l + 2a°o ' 

X> /3S2±SSL 

— V _a o 

A>- A /Ml±2o)£ 
— V -°o 



Jl^, -1/2 < a < 

1 < a < -1/2 
ao < —1 



TABLE II: The real and physically meaningful critical points of interacting model 2 and their behavior. 

In order to make the phase-space analysis more transparent, we elaborate the autonomous system (|26l) numerically 
for the parameter values ao = 0.5 and A = 1.0. The corresponding phase-space trajectories, arising imposing different 
initial conditions, are shown in figure [3l As expected, the system is always attracted by the stable fixed point B of 
table El 



C. Interacting Model 3 

This model is characterized by an interaction term of the form 

where for simplicity we consider n = 1. The autonomous system reads: 

x' = -3x - ^Aj/ 2 + ^x(l - x 2 - y 2 ) - 3(3x(l + x 2 - y 2 ) 



(30) 



(31) 



FIG. 3: Phase-space trajectories for interacting model 2 using Qo = 0.5 and A = 1.0. The stable fixed point is the critical point 
B in table EH with (x c2 , y c2 )={-0Al, 1.08). 



The real and physically meaningful critical points (x c , y c ) of the autonomous system (|3ip are: 



A /A 2 
(x cl = 0, y c i = 0), (x c2 = y C 2 = y 1 + y) 



(»c3, J/c3 = \A - ^c3 - V6Ax c3 /3), (xc4, y c4 = \Jl-x 2 ci - V6Ax c4 /3), 



(32) 



where 



and 



Zc3 



X c 4 



A+y/A 2 -12/3 
2^/3 



A - v/A 2 - 12/3 



2x76/3 



(33) 



(34) 



The matrix S of the linearized perturbation equations is: 



2/c 



A — 3x c 



|Aa; c 



The critical points are presented in table [Till together with their stability behavior, wtot, Q<p and the acceleration 
condition. 



Cr. P. 


X c 


Vc 


Existence 


Stable for 




Wtot 


Acceleration 


A 








all (3, A 


Unstable 








No 


B 


A 

V6 




all (3, A 


all A, P > -1 


1 




all 13, A 


C 


Xc3 


Vc3 


A<^,/?<-l 


Unstable 


-a& + vis 


-Sc3 - y% 


-V^<A< J^,0 < -1 


D 


X C 4 


VcA 


- A /^<A,/3<-l 


Unstable 




2 2 
3-c4 VcA 


-^^<A< v 7 ^, /3<-l 



TABLE III: The real and physically meaningful critical points of interacting model 3 and their behavior. 



In order to present the results more transparently, in fig.2]we depict the range of the 2D parameter space ((3, A) that 
corresponds to the existence of the stable critical point B of table IIII1 Furthermore, we elaborate the autonomous 
system (|3ip numerically, for the parameter values A = 1.0, (3 = 0.5 and A = 1.0, (3 — —2.0 . The phase-space 
trajectories, corresponding to different initial conditions, are shown in figure In these cases, the system is always 
attracted by the stable fixed point B of table Hill as expected. 
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I 







FIG. 4: The range of the 2D parameter space {(5, A) that corresponds to an existing and stable critical point B of table IIIII 
that is in the case of interacting model 3. 




FIG. 5: Phase-space trajectories for interacting model 3 with (5 = 0.5, A = 1.0 (left panel) and f3 = —2.0, A = 1.0 (right panel). 
In both cases the stable fixed point is the critical point B of table Hill with (x C 2, y C 2) = (-0.41, 1.08). 

D. Interacting Model 4 

In the previous subsections we considered interaction terms depending explicitly on the Hubble parameter H, that 
is interaction terms determined by the properties of the whole universe. In the present model we consider the local 
interaction term \2M 



Q = rp„ 



(35) 



where we assume that the coefficient T is a constant. Note that this ansatz, with T > 0, has been used in different 
frameworks to describe the decay of dark matter into radiation (25|, the decay of a curvaton field into radiation (2fj| 
and the decay of superheavy dark matter particles into a quintessence scalar field [27j . 

The complexity of the interaction term (|35[) . comparing to the previous ones, is revealed by the fact that the 
dynamical evolution cannot be elaborated using only the two auxiliary variables (|16p . Indeed, H cannot be eliminated 
from the autonomous equations. In order to achieve that we introduce an additional auxiliary variable as 



Ho 
H ' 



(36) 



where H is a constant. Thus, using (|3"S"|) , the cosmological equations (0), ©, and ([5]) are transformed to the 



9 



following autonomous system: 



-3* + -x 2 - y 2 ) -J^Xy'-^il + x 2 - y 2 ) 



2x 



?/ = -x 2 - y 2 ) - J ^Xxy 
3 

v' = -v(l - x 2 - y 2 ), 
where we have introduced the dimensionless coupling constant 

r 

The 3x3 matrix S of the linearized perturbation equations is: 

-2 i „,2 



2 + 6rr2 V 1 ' x c Vc) 3 



Vc 



Vc 



A — 3a 



2 

3x c v, 



-v^A-3x c + ^ 
-3y c 2 - 

-3y c w c |(i-x2-y c 2 ) 



o 



(37) 



(38) 



In this model there is only one real and physically meaningful critical point presented in table IIV1 together with its 
stability behavior, w to t, &4> and the acceleration condition. 



Cr. P. 


Xc 




Vc 


Existence 


Stable for 




Wtot 


Acceleration 


A 


A 

V6 







all 7, A 


all 7, A 


1 


-1-^ 


all 7, A 



TABLE IV: The real and physically meaningful critical point of interacting model 4 and its behavior. 




-0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 



FIG. 6: Phase-space trajectories for interacting model 4 with A = 0.5. The stable fixed points is the critical point A of table 
HYl with (xc, y c )=(-0.20, 1.02). 

To reveal this behavior more transparently, we solve numerically the autonomous system (|37p for A = 0.5. In order 
to plot the two dimensional phase trajectories for the variables x and y, we project the system onto the v — plane 
and the results are shown in figure O Finally, we mention that for the corresponding model of a canonical field, 
examined in [2l|, using the variable z = Hq/ (Hq + H) the authors found "critical points" for z = 1, which corresponds 
to v = oo in our case. However, as we can see, this point is not a critical point of the system. 



IV. COSMOLOGICAL IMPLICATIONS 



Since we have performed a complete phase-space analysis of various interacting phantom models, we can now 
discuss the corresponding cosmological behavior. A general remark is that this behavior is radically different from the 
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corresponding interacting quintessence models with the same couplings [20, 21]. Furthermore, a common feature of all 
the examined models is that we obtain < — 1 in the whole phase-space, and w to t < — 1 for all the existing and stable 
critical points. Thus, independently of the specific interacting model and of the chosen initial conditions, we will be 
always moving on the same side (below) of the phantom divide. This was expected for phantom cosmology. However, 
we remind that interacting quintessence can conditionally describe the — 1-crossing, in contrast to expectations (l8l . 
EH H3, EH- Finally, we mention that the fact that wtot is not only less than —1/3, as required by the acceleration 
condition, but it is always less than —1 in all stable solutions, leads to H > at all times. Thus, these solutions 
correspond to a super-accelerating universe [301 ] . that is with a permanently increasing H. Ifi/^ooati— >oo then 
we acquire an eternally expanding universe, while if H — > oo at t — > tsn < oo then the universe results to a Big Rip. 
These behaviors are also common in phantom cosmology 0, l3l| . 

A. Interacting Model 1 

In this model, both critical points A and B are stable points, and thus late-time attractors (see table H]). Point 
A, which is stable for a large region of the 2D parameter space (a, A), corresponds to an accelerating universe. In 
addition, since fi^ = 1 it corresponds to a complete dark energy domination. Note that both these features hold 
independently of the a-value if it is positive, which was expected since positive a means positive Q, that is energy 
transfer from dark matter to dark energy. So it leads to a final, complete dark energy domination. In conclusion, 
we see that the stable fixed point A does correspond to an accelerating universe, but cannot solve the coincidence 

problem since it leads to f2dark energy = 1 instead Of fidark energy/fidark matter ~ 0(1). 

On the other hand, the critical point B, for a < — 3 and for A- values that make it stable (note that in this case 
the existence and stability conditions coincide) corresponds to an accelerating universe with < < 1, i.e with 
^dark energy /^dark matter ~ Therefore, for this region of the (a, A) parameter space, the stable fixed point B 

can satisfy the observational requirements, and at late times the universe is attracted by an accelerating solution 
with vacuum and matter densities of the same order. For example, if we take a = —3.1 and A = 0.1 (see left panel 
of fig. [J), then we obtain f2^ = 0.7 and Wtot = —1-03, which is in agreement with current observations within 
95% confidence level [l|. However, as can be seen from figure [TJ the stability region of point B, which produces the 
aforementioned interesting cosmological behavior, is very narrow. Therefore, the late-time attractor A has significantly 
larger probability to be the universe's late-time solution. In other words, the phantom interacting model at hand can 
indeed satisfy the cosmological requirements, but for a very narrow area of the parameter space. 

B. Interacting Model 2 

In this case we observe that the critical points A, C and D are unstable, and therefore they cannot be late-time 
cosmological solutions. The only relevant critical point is B, which is a stable fixed points for the whole (ao,A) 
parameter space. As can be seen from table HIl it corresponds to an accelerating universe with f2^ = 1, i.e., to a 
complete dark-energy domination. Thus, it cannot act as a candidate for solving the coincidence problem. 

C. Interacting Model 3 

In this model we see that the critical points A, C and D are unstable and thus non-relevant. Critical point B is the 
only cosmologically interesting one, since it corresponds to a late-time attractor for the region of the A) parameter 
space shown in figure [H However, although for that region it always gives rise to acceleration, it possesses £1^ = 1. 
Therefore it cannot solve the coincidence problem. 

D. Interacting Model 4 

In this case, the critical point A is the only real and physically meaningful critical point, and it is stable for the 
whole 2D parameter space (7, A). It corresponds to an accelerating universe, completely dominated by dark energy 
(£1$ = 1). Thus this phantom interacting model cannot solve the coincidence problem. 
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V. CONCLUSIONS 

In this work we performed a detailed phase-space analysis of various phantom cosmological models, where the dark 
energy sector interacts with the dark matter one. Our basic goal was to examine whether there exist late-time scaling 
attractors, corresponding to accelerated universe and possessing fidark energy/^dark matter ~ 0(1), thus satisfying the 
basic observational requirements. We investigated four different interaction models, including one with a local, i.e. 
ii-independent, form of the interaction term (interacting model 4). We extracted the critical points, determined 
their stabilities, and calculated the basic cosmological observables, namely the total equation-of-state parameter wtot 
and f^dark energy (attributed to the phantom field) . The key point of solving the coincidence problem is that the 
universe is in the attractor now and it has normal radiation dominated and matter dominated eras before it reached 
the attractors. Note that as long as the interaction term is not too strong, the standard cosmology can be always 
recovered. Once it is in the attractor, the results do not depend on the initial conditions. Thus, one can switch on 
the interaction and consider as initial conditions the end of the known epochs of standard Big Bang cosmology, in 
order to avoid disastrous interference. 

In all the examined models we found that stable late-time solutions do exist, corresponding moreover to an accel- 
erating universe. This feature was expected since phantom cosmology has been constructed in order to always satisfy 
this condition. Indeed, for all the studied models, we did not find any non-accelerating stable solution. However, in 
almost all the cases the late-time solutions correspond to a complete dark energy domination and thus are unable 
to solve the coincidence problem. The only case in which this is possible is in interacting model 1, if we select the 
parameter values from a very narrow region of the 2D parameter space (see figure [TJ . 

In conclusion, we see that the examined models of interacting phantom cosmology can produce acceleration (which 
is "embedded" in phantom cosmology in general) but cannot solve the coincidence problem, unless one imposes a 
form of fine-tuning in the model parameters. This result has been extracted by the negative-kinetic-energy realization 
of phantom, which does not cover the whole class of phantom models, but since it is a qualitative statement it 
should intuitively be robust for general interacting phantom scenarios, too. An alternative direction would be to 
consider more complicated interaction terms, suitably constructed in order to solve the coincidence problem. But 
the interaction term was introduced in phantom cosmology in order to solve the coincidence problem in a simple 
and general way, avoiding the assumptions and fine-tunings of conventional cosmology. Although promising and with 
many advantages, interacting phantom cosmology needs further investigation. 



[1] A.G. Riess et al. [Supernova Search Team Collaboration], Astron. J. 116, 1009 (1998); A. G. Riess et al. [Supernova Search 
Team Collaboration], Astrophys. J. 607, 665 (2004); S. Perlmutter et al. [Supernova Cosmology Project Collaboration], 
Astrophys. J. 517, 565 (1999); D. N. Spergel et al., Astrophys. J. Suppl. 148, 175 (2003); S. W. Allen, et al., Mon. Not. 
Roy. Astron. Soc. 353, 457 (2004). 

[2] S.M. Carroll, V. Duvvuri, M. Trodden, M.S. Turner, Phys. Rev. D 70, 043528 (2004); T. Chiba, Phys. Lett. B 575, 1 
(2003); S. Nojiri, S.D. Odintsov, Phys. Rev. D 68, 123512 (2003); C.G. Shao, R.G. Cai, B. Wang, R.K. Su, Phys. Lett. B 
633, 164 (2006). 

[3] S. Capozziello, Int. J. Mod. Phys. D 11, 483 (2002); S. A. Appleby and R. A. Battye, Phys. Lett. B 654, 7 (2007); S. 

Nojiri, S.D. Odintsov, Int. J. Geom. Meth. Mod. Phys. 4, 115 (2007); A. A. Starobin sky, JETP Lett. 86, 157 (2007); W. 

Hu, I. Sawicki, Phys. Rev. D 76, 064004 (2007); S. Nojiri, S.D. Odintsov, |arXiv:0807 .0685 Irep-th]]. 
[4] G.R. Dvali, G. Gabadadze, M. Porrati, Phys. Lett. B 485, 208 (2000); C. Deffayet, G.R. Dvali, G. Gabadadze, Phys. Rev. 

D 65, 044023 (2002); Y.G. Gong, C.K. Duan, Class. Quantum Grav. 21, 3655 (2004); F. K. Diakonos and E. N. Sari- 

dakis, JCAP 090 2, 030 (2009); Y.G. Gong, C.K. Duan, Mon. Not. Roy. Astron. Soc. 352, 847 (2004); Y.G. Gong, 

|arXiv:0808.1316[ astro-ph]] . 

[5] P. Binetruy, C. Deffayet, D. Langlois, Nucl. Phys. B 565, 269 (2000); R.G. Cai, Y.G. Gong, B. Wang, JCAP 0603, 006 
(2006); Y.G. Gong, A. Wang, Class. Quantum Grav. 23, 3419 (2006); Y.G. Gong, A. Wang, Q. Wu, Phys. Lett. B 663, 
147 (2008); M. R. Setare and E. N. Saridakis, Phys. Lett. B 670, 1 (2008); M. R. Setare and E. N. Saridakis, JCAP 0903, 
002 (2009). 

[6] B. Ratra and P. J. E. Peebles, Phys. Rev. D 37, 3406 (1988); C. Wetterich, Nucl. Phys. B 302, 668 (1988); A. R. Liddle 

and R. J. Scherrer, Phys. Rev. D 59, 023509 (1998); I. Zlatev, L. M. Wang and P. J. Steinhardt, Phys. Rev. Lett. 82, 896 

(1999); Z. K. Guo, N. Ohta and Y. Z. Zhang, Mod. Phys. Lett. A 22, 883 (2007). 
[7] R. R. Caldwell, Phys. Lett. B 545, 23 (2002); R. R. Caldwell, M. Kamionkowski and N. N. Weinberg, Phys. Rev. Lett. 91, 

071301 (2003); S. No jiri and S. D. Odintso v, Phys. Lett. B 562, 147 (2003); V. K. Onemli and R. P. Woo dard, Phys. Rev. 

D 70, 107301 (2004) [arXiv:gr-qc/0406098| ; M. R. Setare, Eur. Phys. J. C 50, 991 (2007); E. N. Saridakis, |arXiv:0811.1333l 

[hep-th]]. 

[8] B. Boisseau, G. Esposito-Farese, D. Polarski and A. A. Starobinsky, Phys. Rev. Lett. 85, 2236 (2000); S. Nojiri, 
S. D. Odintsov and M. Sasaki, Phys. Rev. D 71, 123509 (2005); M. z. Li, B. Feng and X. m. Zhang, JCAP 0512, 



12 



002 (2005); S. Nojiri and S. D. Odintsov, Phys. Rev. D 72, 023003 (2005); S. Sur and S. Das, JCAP 0901, 007 (2009); 
K. Bamba, C. Q. Geng, S. Nojiri and S. D. Odintsov. larXiv:0810.4296l [hep-th]. 



[9] C. Armendariz-Picon, V. F. Mukhanov and P. J. Steinhardt, Phys. Rev. D 63, 103510 (2001) [arXiv:astro^ph/ 0006373 . 
[10] B. Feng, X. L. Wang and X. M. Zhang, Phys. Lett. B 607, 35 (2005); Z. K. Guo, et al, Phys. Lett. B 608, 177 (2005); 

M.-Z Li, B. Feng, X.-M Zhang, JCAP, 0512, 002 (2005); B. Feng, M. Li, Y.-S. Piao and X. Zhang, Phys. Lett. B 634, 101 

(2006); M. R. Setare, Phys. Lett. B 641, 130 (2006); W. Zhao and Y. Zhang, Phys. Rev. D 73, 123509 (2006); M. R. Setare 

and E. N. Saridakis, Phys. Lett. B 671, 331 (2009). 
[11] A. G. Cohen, D. B. Kaplan and A. E. Nelson, Phys. Rev. Lett. 82, 4971 (1999); P. Horava and D. Minic, Phys. Rev. 

Lett. 85, 1610 (2000); S. D. H. Hsu, Phys. Lett. B 594, 13 (2004); M. Li, Phys. Lett. B 603, 1 (2004); Y.G. Gong, Phys. 

Rev. D 70, 064029 (2004); Y.G. Gong and J. Liu, JCAP 0809, 010 (2008); D. Pavon and W. Zimdahl, Phys. Lett. B 

628, 206 (2005); H. Li, Z. K. Guo and Y. Z. Zhang, Int. J. Mod. Phys. D 15, 869 (2006); M. R. Setare, J. Zhang and X. 

Zhang, JCAP 0703, 007 (2007); E. N. Saridakis, Phys. Lett. B 660, 138 (2008); E. N. Saridakis, JCAP 0804, 020 (2008); 

E. N. Saridakis, Phys. Lett. B 661, 335 (2008). 
[12] A. P. Billyard and A. A. Coley, Phys. Rev D 61, 083503 (2000); J. P. Mimoso, A. Nunes and D.Pavon, Phys. Rev. D 

73, 023502 (2006); R. Lazkoz and G. Leon, Phys. Lett. B 638, 303 (2006); T. Gonzalez, G. Leon and I. Quiros, Class. 

Quant. Grav. 23, 3165 (2006); G. R. Farrar and P. J. E. Peebles, Astrophys. J. 604, 1 (2004); B. Wang, Y. G. Gong and 

E. Abdalla, Phys. Lett. B 624, 141 (2005); M. R. Setare, Phys. Lett. B 642, 1 (2006). 
[13] Z. K. Guo, R. G. Cai and Y. Z. Zhang, JCAP 0505, 002 (2005); A. Nunes, J. P. Mimoso and T.C. Charters, Phys. Rev. D 

63, 083506 (2001); D.F. Mota and C. van de Bruck, Astron. Astrophys. 421, 71 (2004); M. Manera and D.F. Mota, Mon. 

Not. Roy. Astron. Soc. 371, 1373 (2006); N.J. Nunes and D.F. Mota, Mon. Not. Roy. Astron. Soc. 368, 751 (2006); J.D. 

Barrow and T. Clifton, Phys. Rev. D 73, 103520 (2006); T. Clifton and J.D. Barrow, Phys. Rev. D 73, 104022 (2006); T. 

Clifton and J.D. Barrow, Phys. Rev. D 75, 043515 (2007); M. Jamil and M.A. Rashid, arXiv: 0802.1144; M. Jamil, arXiv: 

0810.2896. 

[14] P.G. Ferreira, M. Joyce, Phys. Rev. Lett. 79, 4740 (1997); E.J. Copeland, M. Sami, S. Tsujikawa, Int. J. Mod. Phys. D 

15, 1753 (2006); Y.G. Gong, A. Wang, Y.Z. Zhang, Phys. Lett. B 636, 286 (2006). 
[15] E. J. Copeland, A. R. Liddle and D. Wands, Phys. Rev. D 57, 4686 (1998). 

[16] P. Singh, M. Sami and N. Dadhich, Phys. Rev. D 68, 023522 (2003); J. G. Hao and X. Z. Li, Phys. Rev. D 70, 043529 
(2004). 

[17] H. Wei and S. N. Zhang, Phys. Rev. D 76, 063005 (2007); M. R. Setare and E. N. Saridakis, JCAP 0809, 026 (2008); 
M. R. Setare and E. N. Saridakis, Phys. Lett. B 668, 177 (2008); M. R. Setare and E. N. Saridakis, |arXiv:0807.3807l 
[hep-th]]. 

[18] C. Wetterich, Astron. Astrophys. 301, 321 (1995); L. Amendola, Phys. Rev. D 60, 043501 (1999). 
[19] H. Garcia-Compean, G. Garc ia- Jimenez, O. Ob regon and C. Ramirez, JCAP 0807, 016 (2008). 
[20] X. M. Chen and Y.G. Gong, |arXiv:0811.f698"! gr-qc]]. 

[21] C. G. Bohmer, G. Caldera-Cabral, R. Lazkoz and R. Maartens, Phys. Rev. D 78, 023505 (2008). 
[22] L. Amendola, M. Quartin, S. Tsujikawa and I. Waga, Phys. Rev. D 74, 023525 (2006). 

[23] Z. K. Guo and Y. Z. Zhang, Phys. Rev. D 71, 023501 (2005). T. Gonzalez and I. Quiros, Class. Quant. Grav. 25, 175019 
(2008). 

[24] R. Curbelo, T. Gonzalez, G. Leon and I. Quiros, Class. Quant. Gra v. 23, 1585 (2006). 

[25] R. Cen, Astrophys. J. 546, L77 (2001) [arXiv:astro-ph/0005206] ; M. Oguri, K. Takahashi, H. Ohno and K. Kotake, 

Astrophys. J. 597, 645 (2003). 
[26] K. A. Malik, D. Wands and C. Ungarelli, Phys. Rev. D 67, 063516 (2003). 
[27] H. Ziaeepour, Phys. Rev. D 69, 063512 (2004). 

[28] R. J. Scherrer and A. A. Sen, Phys. Rev. D 77, 083515 (2008); R. J. Scherrer and A. A. Sen, Phys. Rev. D 78, 067303 

(2008); M. R. Setare and E. N. Saridakis, Phys. Rev. D 79, 043005 (2009). 
[29] L. P. Chimento, A. S. Jakubi, D. Pavon and W. Zimdahl, Phys. Rev. D 67, 083513 (2003). 

[30] S. Das, P. S. Corasaniti and J. Khoury, Phys. Rev. D 73, 083509 (2006); M. Kaplinghat and A. Rajaraman, Phys. Rev. D 
75, 103504 (2007). 

[31] F. Briscese, E. Elizalde, S. Nojiri and S. D. Odintsov, Phys. Lett. B 646, 105 (2007). 



